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Introduction 

Gauge invariant descriptions of hadrons and their substructures are needed to directly 
link phenomenology to QCD. The need for gauge invariance becomes clear when one 
thinks of jet formation: here a quark and an anti-quark move away from each other and 
turn into two independent jets of hadrons. To obtain such colourless objects, without 
making the quark matter directly interact so that jet formation becomes implausible, one 
needs to surround the quarks with glue, which makes the system locally gauge invariant. 
It is then energetically favoured for this glue to split into two as the quark matter fields 
reach a large enough separation. This is naively described by parton fragmentation [Q] 
or, closer to QCD, modelled by a gluonic string snapping [^. In reality, of course, the 
glue will have a much richer structure. 

It is known |Q that colour charge is only defined on locally gauge invariant states 
and thus all physically meaningful hadronic substructures - coloured quarks in jet for- 
mation, constituent quarks, the form of the glue in glueballs and the pomeron - must 
be described in a gauge invariant way. The appropriate variables for short and long 
distance physics will vary, however, there are far too few methods known with which 
they can be described. In this paper we want in particular to show how to proceed from 
string-based to constituent-type descriptions. We will argue that these latter variables 
provide a better starting point for probing the above structures. 

Non-local states and currents which are supposed to overlap with hadrons are gener- 
ally made gauge invariant with the insertion of 'strings', i.e., the fundamental fields are 
linked by path ordered exponentials. To illustrate this, consider the case of Quantum 
Electrodynamics (QED) where we can make a candidate description of positronium by 
attaching a string that runs from one fermion to the other over a path F as depicted in 
Fig. g. 




Figure 1: Using strings to try to describe e"'"e states. 

This system (which we for simplicity keep in a single time slice) 

^(t, x) exp I -ie J dziAi{t, z) | tp{t, y) , (1) 

is gauge- invariant, however, it has a manifest path-dependence which has no obvious 
physical interpretation. It is, though, possible to remove this path dependence using 



the standard decomposition oi Ai into longitudinal (Af) and gauge-invariant, transverse 
(Af) degrees of freedom: A^ = Af + Af where Af = didjAj/V"^. We so obtain 

V^(x)e'"%^(")e-'"/r'^^'^nM)g-ie%^(j/)^^^^, ^ ^2) 

where V~^ is the standard Green's function for the spatial Laplacian. 

Note that by writing the positronium state in this way we have put all of the arbitrary 
F-dependence into a separately gauge invariant factor and it is natural to drop it. The 
remaining gauge-invariant combination 
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^D{x)tpD{y) ■= ^{x)e^^''''^e *''"^^^V(2/) , (3) 

is a product of two individually gauge invariant terms, ipD, which were first written 
down by Dirac [Q and have been employed by a variety of authors in the intervening 
period (for a review, see p). 

Further strong physical reasons for dropping the path dependent factor in (0) can be 
seen by considering the energy of the (straight line) string state. To find the expectation 
value of the Hamiltonian, H = ^ J d?w{Ef{w) + Bf{w)), for such states, we employ the 
standard equal-time commutators, [Ei{x),Aj{y)] = iSij6{x — y). The magnetic field 
commutes with (|1]) but the electric field does not. We straightforwardly find for the 
potential energy of this state @] 

V{r) = pf\0)r, (4) 

where r = \x — y\ and we have dropped the divergent self-energy contributions (the r 
independent terms). This shows a linear, confining potential with a divergent coefficient 
- in QED! It is clear that 'confinement' here is an artefact due to trapping the electric 
fiux in a string and that the divergence of the coefficient is a result of the string being 
infinitely thin (the _L subscript indicates that this divergent factor originates from the 
directions orthogonal to the path of the string). 

Of course this, although gauge invariant, is not physical. It corresponds to an in- 
finitely excited state and is not stable |^ . A discussion of this instability and an anima- 
tion of the decay of such a state can be found at http : //www . if ae . es/~roy/qed . html. 
We can gain further insight from recalling that Gauss' law tells us that diEi = ep, where 
p is the charge density. If we now decompose the electric field into transverse and longi- 
tudinal components, the contribution of the electric field to the above Hamiltonian can 
be so re-expressed as a sum of two positive semi-definite terms, ^Jd^w{Ef + E^ ). 
The integral over the longitudinal part of the electric field is evidently fixed by Gauss' 
law to be the same for any gauge invariant description of a state with a fixed number of 



charges. However, the transverse component, Ej , is not so fixed. For the above 'stringy' 
state it in fact diverges, while for the variable ipD it is zero which minimises the energy. 
Hence the stringy state is unstable and it decays into a product of two charged fields as 
given in (|^). It is easy to repeat the analysis leading to (^) for these constituents and so 
see that the familiar potential of electrodynamics, — e^/47rr, holds for this state. This 
implies that in any string description of this e^e~ state, where we can factor out the 
exponential over the Af variables, we indeed ought to drop this factor and work with a 
product of two ipD variables. 

The string-based description has, we have seen, many unattractive characteristics 
which make it a poor starting point to analyse positronium. The alternative, constituent 
picture has a much more appealing physical base. Its domain of validity is for heavy 
charges since the charges, -^d, have to be static f^. To approach a more reahstic 
description, one would need to expand in the masses and so build in binding energies. 
However, since this starting point is energetically favoured and has a clear physical 
interpretation it is a more sensible starting point for studying bound states in QED. 

The question now arises, can a similar decomposition be carried out for QCD? The 
string picture, since the matter fields would feel confining interaction at short distances, 
does not include the short distance forces felt by quarks and gluons. However, the 
generalisation of ipD would include the Coulombic type of interaction which dominates 
at short distances. It is not, however, initially obvious that this can be done. The 
vector boson fields are now Lie algebra valued and any T" matrices sandwiched be- 
tween spinors cannot be extracted. Furthermore the decomposition into longitudinal 
and gauge- invariant, transverse degrees of freedom used above does not hold in QCD. 
In the rest of this paper we will show that these obstacles can in fact be overcome at 
any order in perturbation theory and the path dependence can be removed in QCD. In 
the mesonic case we will thus regain the constituent structure which has been previously 
shown 1^ by us to be responsible for asymptotic freedom. 

Non-Abelian Strings 

We now want to recall how to introduce into a non-abelian theory variables which are 
both gauge invariant and do not have a path dependence. Then we will show the relation 
between these degrees of freedom and the string-based systems at lowest order in the 
coupling. 

We define h'^ to be a field dependent element of the gauge group that transforms 
as 

h-^{x)^h-\x)U{x). (5) 

Recalling how fermions transform under gauge transformations, if) -^ U~^ip, it is clear 
that the following fields 

M^) = h-'^ix) , (6) 



describe gauge- invariant quarks. This is the non-abehan analogue oiipD- We call h ^ a. 
dressing, it corresponds here to the gluonic cloud which surrounds a static fermion. 

To make this explicit, we now recall that there is a simple generalisation of (|]) to all 
orders in perturbation theory (see the appendix of |^ for a discussion of how to generate 
it). Essentially we have 

h-' = eM{9Xt + 9\2 + 9\l + ■ ■ ■)T''} , (7) 

where 

X? = ^ , X^2 = r^4 (xiA'j + kd.xDx'i) , etc . (8) 



Here Xi is, up to a trivial colour factor, what we had above in QED while X2 etc are 
minimal extensions whose role is to preserve the gauge transformation property, Eq. ^. 

In a previous letter [§] we have demonstrated that precisely this gluonic configuration 
in h~^ is responsible for antiscreening. We now want to show how we can use these 
descriptions to remove the string dependence in the QCD equivalent of Fig. 0: 

^(a;)Pexp (^g J^ dziAiz)^ i){y) , (9) 

which is gauge invariant but path dependent. At lowest order in perturbation theory in 
QCD we can make the same decomposition of the vector potential, since Af is gauge 
invariant up to terms of order g. We can thus rewrite (|^) as 

^(x)/i(x)Pexp [gjyz,A^{z)^ h-\y)i^{y) + 0(/) , (10) 

where h~^ is the static dressing. We cannot now factor off the exponential due to the 
presence of the Gell-Mann matrices. However, if we write out the last equation more 
explicitly: 
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l + ^^^(y)V(y) + 0(/), (11) 



then, since at this order the terms in the square bracket are gauge invariant, we may 
discard the path dependent term, keeping only the 1, to obtain, at this lowest order, a 
result equivalent to that which we would obtain from multiplying two dressed quarks, 
^/'/j, together. Eq. [TT] is suggestive of how this will generalise to higher orders. We stress 
again that the term which we drop has a singular energy since the flux is concentrated 
onto a line. 

Before moving on to the general proof, let us consider the baryonic case at 0{g). 
Essentially we want to show that the structure of Fig. a can be rewritten in pertur- 
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Figure 2: Two descriptions of a baryon: a) using strings; b) in terms of constituents. 
Note that a) also possesses a dependence on the point where the strings meet and that 
only in b) are there locally gauge invariant substructures with well-defined colour. 



bation theory as a product of three gauge invariant quarks, each then with well-defined 
colour, plus extra additive terms which contain the (singular) path dependence and 
should be dropped. 

The fermions in local baryonic currents^ such as e°'^'^%lj°'ip^ip'^ may be taken to different 
positions in a gauge invariant way by attaching strings as follows 



PexpU/ dziAi{z)\ ^{x) Pexpig dziAi{z)] ipiy) 
Pexpig J dz^Ai{z)] ij{l) 



(12) 



If we decompose the vector potentials once more into transverse and longitudinal com- 
ponents we obtain at 0{g) 



h{w) (l + g r dziAj{z)] h-\x)^{x) 



X 



hiw){l + g dz'.A Uz') h-\y)^iy) 



X 



h{w)il + g dz\Al{z")]h-\l)ij{l) 



+ Oig')e 



2\ nbc 



(13) 



Again discarding the singular path dependent terms, we obtain at this order two types 
of terms: those corresponding to an (antisymmetric) product of three dressed quarks at 
X, y and / and also terms of the form 



arixyii-ii/) 
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e"'" + 2 similar , 



(14) 



^The exact Lorentz structure of the current is unimportant. 



which come from expanding the h{w) factors in (|13[). These terms are not path depen- 
dent, but they appear to have an unattractive dependence upon the arbitrary point, w, 
where the strings join. However, we can show that they vanish; to do this rename the 
dummy indices to extract an overall factor of 



(T"),,e-''- + {T"),be''"' + {T^)cae""' , (15) 

which vanishes due to the tracelessness of the Gell-Mann matrices. We thus see that on 
extracting and discarding the path dependence we may, at lowest order in the coupling, 
write (|T2D as a product of three dressed quarks, iph, each of which has a well-defined 
colour. This system then depends neither on the arbitrary paths nor on their junction 
point. 

It should be clear that we could also attach a string to an individual quark or 
connect field strengths to produce non-local descriptions of glue. We will return to 
these applications below. Now though we wish to show how these arguments can be 
extended to all orders. 

All Orders Strings 

We now want to show how to express path ordered exponentials over the vector potentials 
of QCD in terms of products of path ordered exponentials over the gauge invariant gluons 
and the dressings. This will permit us to remove the contour dependence to any order 
of perturbation theory. 

Recall that a path ordered exponential ^ transforms under a gauge transformation 
as 

Pexp (^gjyziAiz)^ -^ U-\x)Pexp [g j^z^Aiiz)^ U{y) , (16) 

and, since a fermion transforms as ^/^ — > U~^ip, we can produce gauge invariant compos- 
ites by linking quarks together via such strings. 

We now want to demonstrate that we can rewrite this at all orders as 

Pexp [gjyzMz)^ = /i(x)Pexp (^g J' dz,Al{z)^ h-\y) , (17) 

where the external dressing factors are gauge dependent and the central term is a path 
ordered exponential over the gauge invariant field, A^. There are no ordering problems 
for the glue dressed in this fashion. 

It is also clear that (|17D transforms as one would expect of a path ordered exponential: 

/i(x)Pexp (g r dziA'l{z)\ h-\y) -^ U-\x)h{x)P exp (g r dz,A'l{z)] h-\y)U{y) . 

(18) 



To prove ([T7|), recall the formal definition of a path ordered exponential 



N-l 



Pexpfg rdz,Aiz))=\im [] [I + gA.,{z{k)) {z,ik + 1) - z,ik)}] , (19) 

where the end points z(0) and ziN') are fixed at x and y respectively. 

If we consider one element of this product, which we write as 1 + gAi[w){wi — Sj}, 
and take {wi — s^} to be small then we want to show that this term can be written using 
the dressing as 

1 + gA^{w){wi - Si} = h{s) \l + gA'l{w){wi - sjl h-\w) . (20) 



If we can show this for one term, then (|T^ follows by inspection. We write the 1 as 
h{s)h~^{s) and Taylor expand h~^{s) around w. Dropping terms of order {wi — s,}^ we 
so obtain 

h{s) (l + g (h-\s)A,{s)h{s) - -{d,h-\s))h{s)\ {w, - sA h-\w) , (21) 

and using {dih~^{s))h{s) = —h^^{s)dih{s), we obtain (^) and have thus demonstrated 



As well as this formal all orders proof, we have explicitly checked (|T^ at higher 
orders. With this result we will now find it relatively straightforward to extract the 
string dependence from descriptions of hadrons in a gauge invariant way. 

Making Hadrons 

Let us now look at the following in turn: quarks, mesons, baryons and finally glue. 
a) A quark with a string attached to it may be written as 

Pexp [g jyz,A,{z)^ i;{y) = /i(oo)Pexp i^g £dz,Al{z)) h-\y)ij{y) , (22) 

where the string is taken along some path out to infinity. Since for reasons of finite 
energy we require that our fields vanish at spatial infinity, we have h{oo) —>■ 1 and we can 
also factor out the gauge invariant (but contour dependent) path ordered exponential, 
leaving ourselves with a dressed quark, h~^{y)ilj{y). This shows how one can proceed 
from the gauge invariant variables introduced by Mandelstam ||10| to those of Dirac. 



b) A mesonic structure may be made up by linking two fermions by a string 

^(x)Pexp [gjyz,A,{z)^ i,{y) = i,{x)h{x)P exp (^g J^dz^A^iz)^ h-\y)ij{y) , (23) 

7 



where we again have used ([T7|). If we expand the path ordered exponential in g then only 
the lowest term, unity, is independent of the string's path. Since by construction the 
higher order terms are gauge invariant, we can drop them and so obtain the non-local, 
gauge invariant and path independent system, xl!{x)h{x)h~^{y)ip(y). Note that we here 
cannot factor out the F dependence but are free to drop this additive, singular term. 



c) For baryonic structures we recall (0) which, using our result ([T7|) and discarding as 
in the mesonic case all the path dependent terms, can be expressed as 

h{w)h-\x)^{x)y [h{w)h-\y)^{y)f [h{w)h-\z)^{l)Y e^^\ (24) 

which can be rewritten as 

hiw^''' h{wf^' h{wy'' [ij\x)] "' [i'''{y)] ^' [^^/)] '' e'^*" . (25) 

Since both the entire expression and the individual dressed quarks are by construction 
gauge invariant, h{wY°' h{w)^^ h{wY'^ e"'"^ must also be gauge invariant. However, we 
can always choose a gauge such that h{w) = 1, and so hiwY"" h^wY^ h{wY^ e"**^ = 
^aa ^bb ^cc ^abc ^^ g^jj gaugcs. Therefore (|25|) can be written solely in terms of dressed 
quarks as ip'^{x)ipj^{y)ip^{l)e"-^'^. 

d) As far as glue is concerned we can similarly link field strengths with path ordered 
exponentials and take traces. The path dependence can then be extracted just as above. 
Of course we also have for glue the possibility to just directly combine dressed vector 
potentials, A'^A'^^. Such gauge invariant structures have not previously been studied. 
Condensates of such structures, non-local, gauge invariant and of dimension 2 could be 
phenomeno logical significant; we note that many authors have argued for the importance 
of (y4^) condensates - it seems that the stigma of gauge dependence can be removed 
from such objects (but see below). 

Conclusions 

We have seen in QED that string-like descriptions are singular and that the path de- 
pendence can be removed, in a gauge invariant way, to yield the correct lowest energy 
state for static charges. In QCD we have demonstrated that such path dependence can 
be removed at all orders of perturbation theory. This opens up new, gauge invariant 
ways to describe hadrons and their substructure. 

Fig. ^ shows three possible ways of describing a meson in a gauge invariant fashion. 
The approach presented in this paper is closest to Fig. ^b which should describe the 
physics prior to the onset of confinement. (Note though that the fall-off of our variables is 
such that they extend out to spatial infinity.) In particular they include the Coulombic- 
type interaction between quarks. It is thus evident that, at short distances, where 
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Figure 3: Three different descriptions of a meson, a) Strings: suitable at strong coupling; 
b) gauge invariant constituents: describes the short distance glue around current quarks 
and appropriate before confinement dominates; c) hadronic variables: correct at larger 
separations. 

the forces are essentially Coulombic, the descriptions we suggest are to be preferred. 
We note that these fields have been shown to describe the glue responsible for anti- 
screening 0. It can be shown that the constituent picture presented here does not yet 
generate three-body forces in baryons. Such forces may emerge from hybrid pictures 
involving diquarks. 

Infinitesimally thin string-type descriptions are only intuitively expected to be rele- 
vant in the strong coupling limit, where the shortest path is to be used. (We note that 



in weak coupling lattice simulations it is found []TT| in practice that 'smearing' the string 



paths in Wilson loop calculations helps in numerical calculations.) 

Of course, a physical meson, while certainly not well- described by a string, is not 
just two constituents and the favoured glue will be in some sense cigar shaped. The 
emergence of such a picture needs to be further pursued; higher order perturbative 
studies will be useful in this regard, but non-perturbative work will be crucial. We note 
here previous attempts to construct, in simpler systems, gluonic distributions which 
yield a phenomenological potential, of r + 1/r type [|T^] . 

There are many questions which immediately cry out to be addressed. The con- 
struction presented above holds to all orders of perturbation theory, yet we know that 
constituent variables only have a limited phenomenological validity. One of the attrac- 
tive features of our approach is that, although some non-perturbative effects can be 
incorporated into the construction of the dressing, the construction of /i~^ is known to 
break down in the non-perturbative domain ||^. Mapping out its domain of validity is 
a major focus of our programme. An immediate corollary is that the {A'^A'^) conden- 
sate can only have a limited validity. The physical implications of this fact need to be 
studied. 
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We also note that path-ordered exponentials are widely used in the literature to 
create non-local, gauge invariant states. We believe that the path independent approach 
presented here offers an opportunity to test the formalism-dependence of the conclusions 
of such analyses. From what we have seen above, it seems natural to use these variables 
as input for variational methods to find the energetically favoured glue in hadronic 
physics. Finally, the variables used here describe the glue close to current quarks and 
are suitable for short distance hadronic physics. Parton-hadron duality should we feel 
be couched in these terms. As such the route from partons to jets will be modified by 
including dressings and thus these methods have potential implications for jet formation. 

Acknowledgements: We thank Emili Bagan, Robin Horan and Shogo Tanimura for 
discussions. 
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